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We examine the role of intrinsic chain susceptibility anisotropy in magnetic field directed self- 
assembly of a block copolymer using in situ X-ray scattering. Alignment of a lamellar mesophase 
is observed on cooling across the disorder-order transition with the resulting orientational order 
inversely proportional to the cooling rate. We discuss the origin of the susceptibility anisotropy, 

Ay, that drives alignment, and calculate its magnitude using coarse-grained molecular dynamics 
to sample conformations of surface-tethered chains, finding Ay 2 x 10~®. From field-dependent 
scattering data we estimate grains of 1.2 /am are present during alignment. These results demon¬ 
strate that intrinsic anisotropy is sufficient to support strong field-induced mesophase alignment and 
suggest a versatile strategy for field control of orientational order in block copolymers. 

PACS numbers: 83.80.Uv,61.30.Vx,64.75.Yz 


Block copolymers (BCPs) undergo self-assembly that 
results in the formation of periodic structures on meso¬ 
scopic length scales. This self-assembly is driven by 
phase separation of chemically distinct segments, with 
the extent of demixing limited by their physical con¬ 
nectivity. The characteristic length scales therefore are 
largely defined by the size, or molecular weight (MW), 
of the polymer chain. Considerable efforts have been de¬ 
voted to developing methods to reliably direct BCP self- 
assembly, i.e. to align BCP domains, in various device- 
or application-relevant geometries and length scales. 

Under appropriate circumstances, magnetic fields can 
dictate the alignment of BCP mesophases in a highly ef¬ 
ficacious manner. Orientational order develops subject 
to the presence of anisotropic field interactions that are 
sufficiently large to overcome thermal forces. The align¬ 
ment response is a function of the anisotropy in magnetic 
susceptibility. Ay = Xy — X-L? where the parallel direction 
is along the axis of highest rotational symmetry. The free 
energy density of the system, e^, is a function of the angle 
ip between the field and this parallel direction, and of the 
field strength, B, Eq. For an ensemble of anisotropic 
objects, €m can be expressed in terms of the scalar orien¬ 
tational order parameter (P 2 (cos(/:?)) = (^(3cos^ ip — 1)) 
that describes the orientation distribution in the system, 
Eq. The energy difference between orthogonal align¬ 
ments is Acm = —AxB‘^ /2/iq. Strong alignment occurs 


when this magnetostatic energy for structurally coherent 
units, grains, is significant compared to thermal motion, 
that is, for |AP^| = |Ae^|V^ ^ ksT, where Vg = 
is the volume of a grain with a typical dimension It 
is therefore apparent that alignment can occur for suit¬ 
ably large grains, field strengths or magnetic susceptibil¬ 
ity anisotropies. 






(X|| cos^ <^ + X± sin^ ip) 

=^4P(P2(cos(^)) 


( 1 ) 

( 2 ) 


Prior work has relied on liquid crystalline (LC) or crys¬ 
talline assembly of rigid moities integrated with the BCP 
to achieve a sufficiently large Ay to drive alignment at 
reasonable field strengths. Prototypical mesogenic units 
such as cyanobiphenyl species have Ay ~ 10“^ (in SI di¬ 
mensionless volume units) [Il[2]. This has been leveraged 
to drive alignment of several systems |3H6]. Ay becomes 
very small for typical BCPs in the absence of mesogenic 
groups. Shape anisotropy notwithstanding, for a lamel¬ 
lar diblock copolymer with volume fractions (pA and 0^, 
the anisotropy, with respect to the lamellar normal, is 
Ax = -{xa-XbY/[{XA l(t>A) + {XB/(t>B)\ where xa and 
Xb are the isotropic susceptibilities of the blocks. On 
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FIG. 1. (a) Chemical structure of PS-b-P4VP. (b) Tempera¬ 
ture resolved SAXS data (25 ^ 265 ^ 25 °C) under 6 T field 
applied vertically. 2-D diffractograms are shown inset at three 
representative temperatures. To^t=253 °C, d-spacing=9.5 
nm. (c) Schematic illustrating degenerate system alignment 
with lamellar normals perpendicular to the field. 


this basis, for a typical symmetric non-LC BCP such as 
poly(styrene-b-4-vinylpyridine), Ax 0(10-10) [318]. 
This results in a significantly smaller driving force for 
alignment, and as a result magnetic field alignment of 
such ‘coil-coir BCPs has not been observed to date. 

Here, we describe in situ X-ray scattering experiments 
of a coil-coil diblock copolymer subjected to high mag¬ 
netic fields. Cooling across the disorder-order transition 
in the presence of the field induces perpendicular align¬ 
ment of the lamellar normals to the field direction, sug¬ 
gesting the presence of magnetic anisotropy of a non¬ 
trivial magnitude. The existence of this anisotropy is 
rationalized in terms of the intrinsic anisotropy of indi¬ 
vidual Gaussian chains and the non-zero ensemble av¬ 
erage of such anisotropies that is a consequence of the 
organization of the chain junctions along the lamellar in¬ 
terface between the blocks. We use molecular dynam¬ 
ics to simulate trajectories of representative chains and 
thereby estimate the intrinsic anisotropy of the system. 

The system is poly(styrene-b-4-vinyl pyridine) (PS-b- 
P4VP), Figj^, of MW 5.5 kg/mol (K) and PS weight 
fraction, /ps'=0.49, obtained from Polymer Source. All 
data presented here are for this 5.5K material. Consis¬ 
tency checks were performed using a well-purified sec¬ 
ondary sample prepared by living anionic polymerization 
with MW=5.2K and /ps=0.50. SAXS was conducted as 
samples were cooled at prescribed rates from 0.1 to 2 
°C/min across the order-disorder transition (ODT), un¬ 
der fields strengths from 0 to 6 T. Further details re¬ 
garding experimental methods and molecular dynamics 
simulations are available in the Supplemental Material. 

Temperature resolved SAXS data during heating and 
subsequent cooling at 0.3 °C/min under a 6 T field is 
shown in Fig. [l]3. The primary peak at g* =0.066 nm”^ 
corresponds to the lamellar period of 27r/g*=9.5 nm and 
Todt =253 °C. 2-D diffractograms at select temperatures 
show the initial non-aligned state at room temperature, 
the high temperature disordered state where only corre¬ 


lation hole scattering is visible, and finally the aligned 
state produced on cooling. To the best of our ability to 
measure it, ± 0.5 °C, Todt is unaffected by the field ap¬ 
plication. The concentration of intensity along the equa¬ 
torial direction in the diffractogram indicates that the 
lamellar periodicity is perpendicular to the field which 
was applied vertically, i.e. along the meridional direc¬ 
tion. The organization of lamellae with surface normals 
perpendicular to the applied field represents a degenerate 
situation. Fig. [^, indicating that Ax < 0. TEM views 
along and perpendicular to the field. Fig. confirm the 
alignment of the microstructure and its degenerate na¬ 
ture. 



FIG. 2. TEM images and associated fast Fourier transforms 
(FFTs) of PS-b-P4VP aligned at 6 T visualized (a) perpen¬ 
dicular to the field direction and, (b) along the field direction. 

2-D SAXS data during passage across Todt are shown 
in Fig. lit- The evolution of order and alignment are vis¬ 
ible in the gradual disappearance of the diffuse and az- 
imuthally uniform intensity due to correlation hole scat¬ 
tering, and the concurrent emergence of Bragg scattering 
concentrated in arcs centered along the equatorial line. 
The azimuthal dependence of the scattered intensity in 
the Bragg peak, /((^), reflects the orientation distribution 
of the lamellar normals in the plane of the diffractogram. 
The full width at half maximum (FWHM) of I{ip) de¬ 
pends both on the cooling rate and the field strength, 
with the best alignment (lowest FWHM) observed for 
the combination of smallest cooling rate, 0.1 °C/min, and 
largest field, 6 T, Fig. 

The high degree of alignment suggests that the sys¬ 
tem has a magnetic anisotropy well in excess of the ex¬ 
pected Ax ^ 0(10“^^). Experiments conducted with 
the custom synthesized secondary samples recovered sim¬ 
ilar results, dispelling any concerns regarding sample pu¬ 
rity. Lamellar BCPs are known to form anisotropic, ellip¬ 
soidal, grains during nucleation and growth, with faster 
growth along the lamellar normal due to interfacial ten¬ 
sion effects HHO]. Magnetic shape anisotropy for such 
grains is of a trivial magnitude as the suspending fluid is 
a disordered melt of identical composition to the grain. 
Eurther, the shape anisotropy would lead to alignment of 
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FIG. 3. (a) 2-D SAXS data at selected temperatures show¬ 
ing the emergence of aligned lamellae during cooling from the 
isotropic, disordered state at 0.3 °C/min at 6 T. (b) Diffrac- 
tograms for samples under 6 T at 245 °C after cooling across 
Todt at different rates indicated, (c) Dependence of azimuthal 
FWHM on field strength (at 0.3 °C/min) and cooling rate (at 
4 and 6 T), on bottom and top x axes respectively, (d) Diffrac- 
tograms at 245 °C after cooling across Todt at 0.3 °C/min at 
indicated field strengths. 


the long axis of the ellipsoid parallel to the field, whereas 
here the observed alignment is orthogonal. 

The origin of the anisotropy can be understood partly 
by recognizing that real Gaussian chains are anisotropic 
entities as they possess finite end-end distances. Polar¬ 
izability anisotropy has been treated in the context of 
freely jointed statistical segments by Kuhn El yield¬ 
ing Ay = 71—72 = |(<^i ~ ol2)R^ / {R^)id where o^i 
and 0^2 are the polarizabilities parallel and perpendic¬ 
ular to the bond joining statistical segments, R^ is the 
squared end-end distance of the chain, and (i?^)o is the 
unperturbed mean-squared end-end distance. The prob¬ 
lem has also been treated in the context of the worm-like 
chain model by Benoit, Weill, et al [12] [13]. While 
individual chains are anisotropic, the orientation of the 
end-end vectors in a melt are uncorrelated so there is 
no net anisotropy. The situation is different however for 
polymer brushes where one chain end is tethered to an 
impenetrable surface. This geometrically imposes a net 
anisotropy in which the orientation distribution of end- 
end vectors, ((/:?), is peaked along the surface normal 

= 0, and in which the width of the distribution narrows 
with increasing areal density of the chains. The physical 
analogy between surface attachment of chain ends for a 
brush, and chain end (junction) localization at the inter¬ 
material dividing surface of a BCP provides a segue into a 
treatment of chain anisotropy in BCPs. Indeed, the op¬ 
tical anisotropy of tethered chains has been considered 
in detail by Lodge and Fredrickson, for BCP melts HU- 
They note the significance of intrinsic anisotropy rela¬ 
tive to form anisotropy in the birefringence of lamellar 
mesophases and highlight the importance of chain ori¬ 


entation due to correlation in the orientation of end-end 
vectors in the ordered state, relative to chain stretching. 

Optical anisotropies have been investigated experimen¬ 
tally by depolarized Rayleigh scattering [ISH 13 and nu¬ 
merically based on bond polarizability data using the 
rotational isomeric states approaches advanced by Na- 
gai [IH], and by Flory and co-workers m- The statis¬ 
tical segment polarizability anisotropy Aa = a\ — 
however is a quantity associated with a hypothetical 
construct, i.e. the statistical segment of Kuhn’s freely 
jointed model. As such Aa cannot be linked readily to 
chemical structure. It has been measured however for 
some polymers using stress-optical coefficients muin]. 
These data are complemented by studies of electric field 
induced birefringence, the Kerr effect, in polymer so¬ 
lutions Eoi mi- The magnetic analog in the Cotton- 
Mouton coefficient for magnetic birefringence provides 
access to the chain magnetic susceptibility anisotropy as 
Cm = {C^njcB'^) ~ AxAq, with Ax interpreted again 
in the context of the chain statistical element mm- 
Precise measurements of magnetic birefringence and seg¬ 
mental susceptibility anisotropy of flexible chains in the 
melt are challenging, and there is little reported data 
that can be used in the present context. Additionally, 
attempts to link statistical segment anisotropy to that of 
the physical monomer lead to mixed results, particularly 
for polymers such as PS due to phenyl ring rotation m- 

In lieu of robust experimental data, we use molecu¬ 
lar dynamics to estimate Ax for PS chains tethered at 
an impenetrable surface. We investigate the orientational 
order of backbone bonds and the side bonds to the phenyl 
ring as a function of degree of polymerization, A/^, and 
the chain areal density at the surface, a. Fig. The 
height of the brush scales as H r\j as expected 

for chains in the melt [24|. The orientational order with 
respect to the surface normal, the y-axis in Fig. [^, for 
backbone and side bonds to the phenyl ring, (Pi (cos (^)) 
and (PI (cos (^)) respectively, increase slowly, but approx¬ 
imately linearly with a. There is no statistically signifi¬ 
cant MW dependence, highlighting the lack of persistence 
of the chains. (P|) and (PI) differ in sign, reflecting the 
tendency of the backbone bonds therefore ^ to lie par¬ 
allel to the surface normal, and the side bonds to be 
orthogonal to the backbone. 

The areal density of chains at the block interface is 
a = pLqNa/MW where Lq is the lamellar period. Here, 
a ^ 1 nm“^ which corresponds to a = 0.033 in simu¬ 
lation units. The orientational order of the phenyl ring 
defined by the ring normal, (Pf (cos 99 )), is also calculated 
(Supplemental Material). At this density, (Pi(cos(^ 9 )) 
0.07 and (P^icosp)) ^ 0.04. We estimate Ax for 
the brush based on the volume fraction, 0, weighted 
anisotropies of the alkane backbone (Ax^) [25] and that 
of the phenyl ring (Ax^) [H]. Our reference direction 
is the surface normal. For the PS chains then. Ax = 
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FIG. 4. (a) Schematic showing a single highlighted 

polystyrene chain tethered at an impenetrable surface. The 
phenyl rings can rotate around the side-bond to the backbone, 
(b) Simulated backbone, {P 2 ), and side chain, {P 2 ), as func¬ 
tions of areal density a for varying Np. (c) Simulated brush 
height, H{a), for varying Np. 


{P 2 {cos ip)) Ax^ + (cos (^)) « —1.6 x 10“^. 

We expect a similar contribution from the P4VP block 
given the near identical susceptibilities of pyridine and 
benzene, and so for the system overall Ax « —1.6 x 10“^. 
(Additional details in Supplemental Material). 

We use the estimated Ax to consider the alignment 
data of Fig. The azimuthal intensity dependence I{(p) 
reflects the probability of observing lamellar normals at 
a given angle cp with respect to the field direction. This 
probability p{ip,B), is governed by a Boltzmann factor 
incorporating the angle-dependent magnetostatic energy, 
Erp{ip^B) = — (P^/2/io)AxC^ cos^ (/9, for a characteristic 
grain size, Eq. The orientational order parame¬ 
ter resulting from this probability distribution may be 
evaluated by integration, Eq. The orientational order 
parameters are obtained as a function of field strength 



FIG. 5. Field-dependent orientational order parameters (Eq. 
0 for various grain sizes, <f, for Ax = —1.6 x at 253 

Black circles show experimentally derived values for sam¬ 
ples cooled 0.3 °G/min. Fitting Eq. (red circle trace) to 
experimental data yields ^ 1.2 /xm. 


from Gaussian fits of the SAXS intensity profiies. We 
estimate the representative grain size by fitting the ex- 
perimentaiiy determined (P 2 ) to those caicuiated using 
Eq. yieiding ^ ^ 1.2 /xm. This corresponds to a fieid 
interaction AE^ > 10^ kT. The resuits are shown in Eig. 

with caicuiated fieid dependent P 2 for different grain 
sizes. It is important to note that the estimated grain size 
oniy provides the characteristic dimensions of a struc- 
turai unit which wouid, at steady state, reproduce the 
orientation distribution measured in the aiigned sampies 
at room temperature. Kinetic effects reiated to changes 
in temperature dependent nucieation and growth rates, 
as weii as any Ostwaid ripening, during the cooling ramp 
do not factor into this analysis. 


p{ip,B) 


(P 2 (cos (/:?)) 


e simp dip 

sin ip dip 

fo I P ~ sin ip dip 

simp dip 


The grain size determined by fitting is in modest agree¬ 
ment with an estimated size of ~ 3 /xm from sampling 
TEM images of field aligned samples. To provide bet¬ 
ter quantification we independently measured the grain 
size using a recently developed ‘variance scattering’ tech¬ 
nique [26]. This approach statistically analyses the az¬ 
imuthal intensity variation of a primary Bragg peak, from 
which the number of independent scatterers, and thus the 
grain size, can be quantified. On this basis we estimate 
^ ^ 5 — 6 /xm at room temperature (further details in 
Supplemental Material). The fact that TEM and scat¬ 
tering derived grain sizes exceed the estimate based on 
the orientation distribution reflects the occurrence of mi¬ 
crostructure coarsening without accompanying improve¬ 
ment in the alignment under the field. We should there¬ 
fore view the estimated ^ as an upper bound for the char- 

















5 


acteristic grain size that pertained during the period in 
which the system was responsive to the field. 

BCP nucleation and growth kinetics have been ex¬ 
plored experimentally and theoretically with some suc¬ 
cess The slow dynamics of high MW BCPs 

has largely restricted studies to unentangled melts and 
concentrated solutions as a matter of convenience. Grain 
sizes of ^ 0.5 — 5 /im are not uncommon. It is therefore 
reasonable to expect that sufficiently large grains can be 
prepared in a variety of systems to permit control of ori¬ 
entation using magnetic fields. The slow dynamics of en¬ 
tangled melts and the inaccessibility of the ODT in high 
MW systems suggest that this strategy would be limited 
to melts below entanglement MW and concentrated solu¬ 
tions. It is clear that the aromatic nature of PS-b-P4VP 
contributes to a markedly larger Ay than one would en¬ 
counter for other common but non-aromatic polymers 
such as polyethylene oxide or polymethylmethacrylate. 
The fact however that Em ^ means that an order of 
magnitude decrease in Ay can be compensated by an in¬ 
crease in ^ by a factor of just over 2, for the purposes of 
alignment. 
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SUPPLEMENTAL MATERIAL 


SAMPLE PREPARATION 

Poly(styrene-b-4-vinylpyridine) (PS-b-P4VP) with 
molecular weight MW=5.5 kg/mol (2.7K PS/2.8K 
P4VP) and polydispersity index D^ 1.2 was obtained 
from Polymer Source and used as received. A secondary 
sample with MW=5.2K (2.6K PS/2.6K P4VP) was 
prepared by living anionic polymerization. Material was 
pressed into an aluminum sample stage supported by 
thin Kapton windows at 160° C to form cylindrical discs 
a few mm thick, and 1 cm in diameter. Small-angle 
x-ray scattering (SAXS) was performed in vacuum with 
a pinhole-collimated instrument (Rigaku SMAX3000) 
using Gu Ka radiation (A = 1.542 A). The beam has a 
1 mm diameter at the sample plane and the accessible 
range of scattering vectors is 0.02 to 0.35 A“^. SAXS 
data were calibrated using a silver behenate standard. 
2-D patterns were integrated azimuthally to provide 
1-D representations of I(q), where q = (47r/6>) sin6> 
and 20 is the scattering angle. Temperature-resolved 
measurements were performed with a custom made 
hot-stage with an associated Omega temperature con- 
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troller. Samples were loaded onto the temperature 
controlled aluminum stage within a superconducting 
magnet (American Magnetics Inc.) with a flux density 
between 0 and 6 T. Samples were subjected to varying 
cooling rates and field strengths, typically by heating 
to 265°C and then cooling to room temperature in the 
presence of the field. For temperature resolved SAXS 
measurements, data were acquired in a quasi-isothermal 
experiment in which the sample was held at constant 
temperature for several minutes at closely spaced tem¬ 
perature points. For TEM, ^50 nm thin sections were 
prepared using a diamond knife mounted on a Leica EM 
UC7 ultramicrotome. Sections were collected on TEM 
grids after flotation onto water and stained inside an 
iodine chamber for 1 hour prior to imaging. Sectioned 
samples were then visualized by EEI Tecnai Osiris TEM 
with an accelerating voltage of 200 kV. 



FIG. 6. TEM image of PS-b-P4VP (2.7k/2.8k) prepared as 
described above but cooled in the absence of the field at 0.3 
°C/min. Here, we show the lamellar microstructure (9.5 nm 
periodicity) of a representative single grain, >2 /xm x 2 /xm. 
Grains of this size and larger are present throughout the sam¬ 
ple whether or not a magnetic field is applied. 


ODT DETERMINATION 

To calculate the order disorder temperature of PS-b- 
P4VP (2.7k/2.8k) we plot inverse intensity of the primary 
scattering peak, g* = 0.066 nm“^ vs. inverse tempera¬ 
ture, in K. Todt is extracted as the value of temperature 
at = 0 by extrapolating a linear fit of the inverse 
intensity drop during ordering. 


Pitted data: l/I = (-18445.795/T) + 35.066 

^ Todt = 252.9°C 



1/Temperature (K’"') 


FIG. 7. Plot of inverse scattered intensity vs. inverse tem¬ 
perature used for Todt determination. 


GRAIN SIZE DETERMINATION BY VARIANCE 
ANALYSIS 

For this measurement and calculation, PS-b-P4VP 
(2.7k/2.8k) was cooled 0.3°C/min in the absence of the 
field in the same set up as described above. Kapton win¬ 
dows were gently removed, and the sample was mounted 
on a tripod polisher using a hot glue (T^ = 80°C) and 
wet-polished using a descending sequence of diamond lap¬ 
ping films on a rotary disc polisher. Before polishing the 
other side of each sample to its final thickness (100-200 
/xm), the specimens were glued onto a 25 /xm thick Kap¬ 
ton foil. The procedure yielded SAXS specimens with 
smooth and parallel sides of known thickness (t) deter¬ 
mined by a micrometer down to a ± 2/xm accuracy. The 
actual thickness was recorded for each sample to be fur¬ 
ther used in calculation of the scattering volume. 

Small angle scattering experiments were performed on 
a Bruker Nanostar system fitted with a 2048 by 2048 pixel 
Vantex-2000 area detector positioned ^ 100 cm from the 
sample. The system was configured to high resolution 
mode using 300 fim scatterless Ge pinholes. The beam 
profile was assessed by measuring transmitted intensity 
during two orthogonal scans of a metal edge located at 
the sample position across the footprint of the beam. The 
vertical and horizontal profiles of the beam were recon¬ 
structed by differentiating the recorded data and fitting 
them to a Gaussian function. The horizontal and vertical 
full width at half maximum (FWHM) parameters were 
used as a measure of the beam width and height, respec¬ 
tively (both slightly larger than 300/xm). The scattering 
volume is calculated as V = FWHMh x FWHMy x t. 
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SAXS data was collected under vacuum, in transmission 
mode, with each frame collected from an independent 
spot for 1 h. Silver behenate was used as the q-space cal- 
ibrant. After the measurements, the patterns were spa- 
cially unwarped using vendor software, and exported for 
further analysis using custom-written Python routines. 

Two-dimensional diffractograms were azimuthally inte¬ 
grated to yield I vs.q curves. The intensity profile of the 
first order 001 {q^) peak was fitted to a Gaussian function 
to measure its position q, and characterize breadth in ra¬ 
dial direction, <jg. From this we computed the expected 
peak width in the azimuthal direction y (i.e., along the 
arc of the scattering ring), The intensity along the 
001 ring in the angular (y) direction was extracted by 
integrating the full peak width (in the q direction), and 
subtracting the local background scattering just outside 
the peak. This I vs.y data was smoothed using a Gaus¬ 
sian function with a width, This averages out noise 
in the signal (e.g., shot-noise due to counting statistics) 
without averaging out the intensity variation arising from 
the grain structure. We then quantitatively convert this 
intensity variation into an estimate for the numbder of 
independant grains within the scattering volume using a 
previously published method [26]. 


Briefly, this method involved computing the relative 
standard deviation of I vs y, which can be converted into 
an estimate of the number of grains (Ng) using knowledge 
of the peak width (cr^), and the peak multiplicity (m=2 
for lamellar peaks). This analysis implicitly assumes 
an isotropic distribution of grain orientations, with rela¬ 
tively well-defined grain boundaries. The average grain 
size can then simply be computed from the known scat¬ 
tering volume (V) as ^ = (317/47rA^^)^/^, which assumes 
roughly spherical grains. We minimize the error of this 
estimate by combining the results from different sample 
thicknesses. 

MOLECULAR DYNAMICS SIMULATIONS 

We perform molecular dynamics (MD) simulations 
on an array of A = 5x5 = 25 coarse grained 
polystyrene chains with degree of polymerization Np = 
10,15, 20, 25, 30 sitting on a square lattice of spacing, d. 
Each monomer (styrene) has n = S coarse-grained atoms 
including two sp^ backbone carbons and six sp^ phenyl 
ring carbons. The total potential energy as a function of 
atom positions (ri,r2, ...rN) is defined as: 


f/(ri,r 2 ,...rN) = ^ULj{rij) + E '^spring T E '^an5'Ze(l*i9 U ’ ^k) T ^ ^ '^torsion (l*i 5 U ’ I*l) T ^ ^ "^wall (^^i) 


i<j 


<ij> 


<ijk> 


<ijkl> 


( 5 ) 


The first term summing over all pairs of atoms is the 
Lennard-Jones (LJ) potential. 


ULj{rij) = 4e 




(6) 


where e and <Jo sets the unit of energy and length of the 
system respectively. 

The second term, summing over all pairs of atoms that 
are connected by chemical bonds, is the harmonic spring 
potential. 


^ spring 2^^s{'kij ^o) : 


( 7 ) 


where the equilibrium bond length is vq = cto, and the 
bond stiffness kg = 9000 and 18000 e/crg for single and 
double carbon-carbon bonds, respectively. 

The third term, summing over every three atoms that 
form a bond angle 6>, is 


Uangle{ri, Tj, Fk) = Uangle(0) = ^ka{0 - 9of, (8) 

where the equilibrium bond angle = 109.5° for sp^ car¬ 
bons and 120° for sp^ carbons. The bond angle stiffness 
is kg = 400 e/rad^ for sp^ carbons and 18000 e/rad^ for 
sp‘^ carbons. This potential keeps the phenyl ring rigid 
and flat. 

The fourth term, summing over every four consecutive 
backbone atoms or carbon atoms connected by sp^-sp^ 
side chains that form a dihedral torsion angle 0, is 


'^torsion(^^i: U ’ ^0 — '^torsion{4^^ — C 3 COS (0 0o) CiCOs(^(f) 0o) T C 3 Ci 


(9) 


where ci,C3 = 6,10 e and the phase shift = 180° figuration, 
such that the four atoms tend to stay in the trans con- 








The last term describes the purely repulsive LJ wall 
at ^ = 0 where polymer chains are attached. In the 
simulation, the first atom of each chain is pinned at ^ = 0. 
Periodic boundary conditions are applied in the x and z 
directions. 

All simulations are run at the reduced temperature 
T = 3.0 e/kB controlled by the Gaussian constraint ther¬ 
mostat, where ks is Boltzmann’s constant. Polymer 
chains are initialized as straight lines with side chains 
perpendicular to the backbone. We first equilibrate the 
system for 10^ time steps until the energy fluctuations be¬ 
come stable, and then sample thermodynamic quantities 
over 10^ time steps with a time step dt = 0.001 i/mcr^/e. 

We measure the backbone side chain P^, and 
phenyl ring PJ bond order parameters with respect to 
the director normal to the surface, where Ok is the an¬ 
gle between the bond or ring normal, and the director 

y- 



There are N^, = {2Np — 1)/N backbone vectors. Ns = 
NpN side chain vectors, and = NpN phenyl ring nor¬ 


mal vectors. For comparison, we also measure PJ with 
respect to the other two axes x and z parallel to the 
surface. 

The average height H of polymer brushes as a function 
of their surface density a = l/d^ is measured for lattice 
spacing d = 4.0,4.5, 5.0, 5.5, 6.0, 6.5, 7.Ocro- H is defined 
as the average of the maximum height coordinates of in¬ 
dividual chains. The scaling exponent am.H r\j {i/d^r 
is extracted as | as expected for chains in the melt. 

ESTIMATION OF SUSCEPTIBILITY 
ANISOTROPY, Ay 

As expected, there is a geometrical relationship be¬ 
tween the orientational order of the side chain, and 

the orientational order of the phenyl ring, (Pf), Figs. 
and The calculated ratio of the two has a magnitude 
of ~ 1.8 which compares well with the expected value, 2. 
We use the simulated orientational order parameters of 
the backbone and phenyl ring, {P 2 ) and (Pf), to esti¬ 
mate the overall anisotropy. Ay, present in the system. 
Although the anisotropy of benzene is commonly con¬ 
sidered to be positive as it is viewed intuitively with re¬ 
spect to the plane of the ring, here we utilize the formal 
specification which is with respect to the axis of highest 
rotational symmetry, the ring normal, and thus Ay^ < 0. 


Ay = (j)^ {P2 (cosO)) Ax^ + (Pf (cos6>))Ay^ « —1.6 x 10 ^ (11) 


Anisotropy of alkane backbone, Ay^ —5 x 
Anisotropy of phenyl ring, Ay^ —5 x 10“^ 
Orientational order of backbone, (P|(cos^)) ~ 0.07 
Orientational order of phenyl ring, (P^icosO)) ^ 0.04 


Phenyl ring volume fraction (estimated as weight fraction), 

0p (MPhenyl/yfstyrene^ _ 

Backbone volume fraction (estimated as weight fraction), 
(j)b ^ 1 — (j)p ^ 0.25 
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FIG. 8. phenyl ring order parameter, as calculated in 

MD simulation, plotted against areal density of chains, a. 



Areal Density, a 

FIG. 9. Ratio of side chain order parameter, (P|) to phenyl 
ring order parameter, (Pf), plotted against areal density of 
chains, a. 




